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Abstract 

J> ; The bilocal electric field strength correlators in Abelian-projected SU (2)- and SU (3)-theories are 
derived with accounting for the contributions to these quantities brought about by the virtual 
■ vortex loops, built out of the dual Nielsen-Olesen strings. Owing to the screening of magnetic 
charge of the dual vector bosons in the gas of such loops, these bosons become heavier, which 
leads to the respective change of the correlation length of the vacuum in the models under study. 
Besides that, it turns out that such a screening leads also to the appearance of the long-range 
^ ■ contribution to one of the two coefficient functions, which parametrize the bilocal correlator within 
the Stochastic Vacuum Model. Finally, the modifications of the propagators of the dual vector 

Ph. bosons, inspired by the correlation effects in the gas of vortex loops, are also discussed. 
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1 Introduction 



Stochastic Vacuum Model (SVM) p] is nowadays commonly recognized to be one of the most 
promising nonperturbative approaches to QCD (see Ref. @ for recent reviews). Within the 
so-called bilocal or Gaussian approximation, well confirmed by the existing lattice data |J, 
this model is fully described by the irreducible bilocal gauge-invariant field strength correlator 
(cumulant), ((F fiu (x)^(x, x')F\ p (x')<&(x' , x))). Here, = d p A v — d v A p — ig[A^,A u ] stands for 

the Yang-Mills field strength tensor, = 4-Pexp ig J A p (u)du p is a parallel transporter 

factor along the straight-line path, and ((OO'}} = (OO'} — (O) (O'} with the average defined w.r.t. 
the Euclidean Yang-Mills action. It is further convenient to parametrize the bilocal cumulant by 
the two coefficient functions |TJ, Q as follows: 

— ((F^(x)$(x, x')F\ p {x')${x' , x))) = 1 Nc xn c {{5^\5 vp - 5p P 5 v \)D [{x - x') 2 ) + 



+- 



d x {{x- x') x 5 up - (x - x') p 5 vX ) + d£((x- a/), 



Di ((x-x') 2 ) 



After that, setting for the nonperturbative parts of the D- and Di-function various Ansatze, one 
can employ SVM for precise calculations of the processes of high-energy scattering || or test 
these Ansatze in the lattice experiments || However from the pure field-theoretical point of 
view, it remains a great challenge to derive the coefficient functions analytically. Unfortunately, 
in this way no big progress has up to now been achieved in the QCD itself. There have rather 
been derived some relations between cumulants of various orders [0], which might be useful only 
in testing the IR asymptotic behaviours of the coefficient functions. 

Contrary to that, more progress has recently been achieved in a derivation of the bilocal 
cumulant in some models allowing for an analytical description of confinement. Those include 
Abelian-projected (AP) SU{2)- 0] and SU(3) |§ theories, as well as compact QED ||. The 
bilocal field strength cumulant in these theories has been studied in Refs. 
(see Ref. 



10, P, pj, respectively 



131 for a review) 



The present paper also follows this line of investigations and is devoted 
to the improvement of calculations of the bilocal cumulant in AP theories. This improvement is 



based on the well known fact [14] that in the case of zero temperature under study, Abrikosov 



vortices [IS| in the Ginzburg-Landau theory (dual Nielsen-Olesen strings |]TB| in our 4D-case) 
form bound states, built out of a vortex and an antivortex, which are usually referred to as vortex 
dipoles (vortex loops in 4D). Such vortex loops are short living (virtual) objects, whose typical 
sizes are much smaller than the typical distances between them. This means that similarly to 
monopoles in 3D compact QED, vortex loops form a dilute gas. The summation over the grand 
canonical ensemble of vortex loops in such a dilute gas approximation was performed in Ref. |T7| 
for the case of the usual Abelian Higgs model (dual AP 5'[/(2)-gluodynamics) in 3D and 4D 
and then extended to the case of the 4D AP 5'[/(3)-gluodynamics in Ref. ]18[ []. On the other 
hand, in all the investigations of the bilocal field strength cumulants in AP theories, performed 
in Refs. JTO, O, the contribution of vortex loops to the partition function, and consequently to 



the cumulants themselves, was disregarded. As it will be demonstrated in the next Section, this 



^^Note that in the case of the 2D Abelian Higgs model, such a summation has for the first time been performed 
in Ref. 111. 
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approximation is really valid, since it is equivalent to another one, which states that the typical 
sizes of vortex loops are negligibly small. However, such a neglection of the contribution of vortex 
loops makes the calculations of the field strength correlators, performed in the above mentioned 
papers, essentially classical. The improvement of these calculations, presented in this paper, is 
just based on accounting for the correlations in the gas of vortex loops. Clearly, such correlations 
are of the quantum origin, as well as the grand canonical ensemble of virtual vortex loops itself. 
Besides that, we shall also evaluate the contributions of vortex loops to the propagators of the dual 
vector bosons and discuss the so-emerging modifications of the respective classical expressions. 

The paper is organized as follows. In the next Section, we shall firstly review the main aspects 
of a derivation of AP theories, necessary for the future purposes. Secondly, we shall review 
the main results of a calculation of electric field strength correlators in the approximation when 
the contribution of vortex loops to these quantities is disregarded. In Section 3, after a brief 
review of the properties of the grand canonical ensemble of vortex loops, we shall evaluate the 
contribution of these objects to the field strength correlators. In Section 4, the same will be done 
for the propagators of the dual vector bosons. The main results of the paper are summarized in 
Conclusion. In the Appendix, some technical details of the calculation of a certain typical integral 
from the main text are outlined. 



Electric Field Strength Correlators in the Absence of 
Vortex Loops Revisited 



2.1 The Models 

To derive from the gluodynamics Lagrangian the IR effective SU(2)- and SU(3) AP theories, 
based on the assumption of condensation of Cooper pairs of AP monopoles, one usually employs 
the so-called Abelian dominance hypothesis |[20|| . It states that the off-diagonal (in the sense of the 
Cartan decomposition) fields can be disregarded, since after the Abelian projection those can be 
shown to become very massive and therefore irrelevant to the IR region, where confinement holds. 
Let us start our analysis with the SU(2)-theoiy. Then, the action describing the rest, diagonal, 
fields and AP monopoles reads 



S, 



eft'. 



t 



J d x (^f^ u 



(2) 



Here, a. 



= A 3 



fnv = d^a u — d u a^, and the monopole field strength tensor obeys Bianchi 



identities modified by monopoles, d^f™ = ^£^v\ p d^f™ p = j™. The monopole currents j™'s should 
eventually be averaged over in the sense, which will be specified below. 

To proceed with the investigation of the monopole ensemble, it is reasonable to cast the theory 
under study to the dual form. This yields the following expression for the partition function: 



Z 



exp 



- d A x 



4 v 



(3) 



where B„ is the magnetic vector-potential dual to the electric one, a 



// 1 U1 "'6" VI " 1V ' "•" |/V I l lil 1( VI VLIKM IV I1U V 1V\ 11.1A VMIV. U /t , 311(1 FflU '.'// •» ' I ■< V,/J~>f, . 

Once the j™-dependence of the action became explicit, it is now possible to set up the properties 
of the monopole ensemble. To describe the condensation of monopole Cooper pairs, it is first 
necessary to specify as the collective current of N of those: 



d»B v 



d u B u 



3 



N 

2g m J2fdx«(s)5(x-x n {s)). 

n=l 



Here, the world line of the n-th Cooper pair is parametrized by the vector x™(s), and g m is the 
magnetic coupling constant, related to the QCD coupling constant g via the topological quanti- 
zation condition gg m = Ann with n being an integer. In what follows, we shall for concreteness 
restrict ourselves to the monopoles possessing the minimal charge, i.e. set n — 1, although the 
generalization to an arbitrary n is straightforward. Secondly, it is necessary to set for the measure 
(. . .)j m the following expression f2lf : 



exp [i d xB^f 



N=l 



N\ 



N 

n 

n=l 



dSr 



i\rj 2 s n 



Vu(s' 



u(0)=u(s n ) 




L ■ 2/ / 



+ 2ig m u /1 (s 



BM4)) " 4A E Jds'tj dslS [«(afl - . (4) 



!,fc=l 



Here, the vector u^(s' n ) parametrizes the same contour as the vector Clearly, the world-line 

action standing in the exponent on the R.H.S. of Eq. contains besides the usual free part also 
the term responsible for the short-range repulsion of the trajectories of Cooper pairs. Equation ([|) 
can further be rewritten as an integral over the dual Higgs field, describing magnetic Cooper pairs, 
as follows: 



exp [i / d xB^f 



exp 



d x 



-\D,M 2 



x(m 2 



(5) 



where = <9 M — 2ig m B fl is the covariant derivative f\ Finally, substituting Eq. (|5]) into Eq. (^), 
we arrive at the following IR effective AP theory of the >S[/(2)-gluodynamics: 

Z = J |$| V |$| VOVB^ exp |- J d 4 x + 1 \D^\ 2 + A (|$| 2 - V 2 ) 2 J , (6) 

where = |$(x)| e ie ^ x ^. Clearly, as soon as we have disregarded the off-diagonal degrees of 

freedom and demanded the condensation of monopole Cooper pairs, this theory is nothing else, 
but just the dual Abelian Higgs model. 

Analogous considerations can be applied to the S'[/(3)-gluodynamics. The only difference 
is that since the 5 , f/(3)-group has two diagonal generators, the resulting AP theory will also be 
[C/ (l)] 2 magnetically gauge-invariant. Within the Abelian dominance hypothesis, the initial action 
reads 



5*. 



off. 



o f m 

A" 1 fJ.u 



d 4 x (f^ + f t 



(7) 



and after the dualization we have for the partition function (cf. Eq. 



2 A seeming divergency at large proper times produced in Eq. (Q) by the factor e AXri s ™ is actually apparent, 
since the last term in the exponent on the R.H.S. of this equation yields the desired damping. 
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Z = Q PB^expj- J 



d x 



Here, F^ = d^ v — c^B^ is the field strength tensor of magnetic field B M , which is dual to the 
field = (^A^,A^j, and j™ = d^f™. A minor nontriviality, one meets further w.r.t. the simplest 
SU(2)-case, is the necessity to take into account the fact that monopole charges are distributed 
over the lattice defined by the root vectors, which have the form 



ei = (1,0) , e 2 



'_1 _V3\ (_1 V3\ 

,2' 2 J ' 63 ^ 2' 2 J ' 
These vectors naturally appear within the Cartan decomposition of the original set of gluonic 
fields as the structural constants in the commutation relations between the diagonal and so-called 
step (raising and lowering) operators. The collective current of iV monopole Cooper pairs then 
reads 



N 3 



29m EEeJ dx<p n (a)6 (x - x^ n (s] 
— i — i 



n=l a=l 

As far as the average over the currents is concerned, it has the form 



(9) 



exp yi J e^xBjJ* 



x exp 



j 

N s } 



a=l I N=l JV ' 



N rds 

n / 



n=l 



-<> 1A ^ / /)//"•(.<) 

u W(0)=«( a )(s n ) 



f / 1 2 

E / ^ hi (« (a) (^)) +2^ m ^( S ;)e a B, 
'=io V 



i/. 



J,Jfc=l 



-4AE jds'.jdslS [u^is^-u^isl) 
1 




3 



| D$ a P$: exp J - 1 rf 4 x p [- - 2«0 m e a B„) $ a | 2 + A (|$ e 



(10) 



where the vector uff(s' n ) parametrizes the same contour as the vector x^ n (s). Finally, it is worth 
noting that since monopoles are distributed over the root lattice, whose vectors are related to each 

3 

other by the condition J2 e a = 0, the dual Higgs fields $ a 's are also not completely independent of 

3 

each other. In Ref. H, it was argued that the relevant constraint for these fields reads 9 a = 0. 

a=l 

Taking this into account we arrive at the following partition function describing an effective [t/ (l)] 2 
magnetically gauge- invariant AP theory of the S'[/(3)-gluodynamics |§: 



Z = f \$ a \V\$ a \V0 a VBr5 (E e a) 

\a=l I 



x exp | — 
where <3> a = I $ a I e l6a . 



d x 



4 w 



E 

a=l 



- \{d„ - 2ig m e a B ll ) ® a \ + A (\<f> a \ A - rf 



(11) 
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2.2 Bilocal Electric Field Strength Correlators 
2.2.1 ST/(2)-case 

In order to investigate bilocal cumulants of electric field strengths in the models (§) and ([0]), it 
is necessary to extend them by external electrically charged test particles (i.e. particles, charged 
w.r.t. the Cartan subgroup of the original £77(2)- or S77(3)-group). For brevity, we shall call these 
particles simply "quarks" . In the SU (2)-case, such an extension can be performed by adding to 

the action (0) the term i J d^xa^j* with j e Ax) = g § dx^(s)5(x — x(s)) standing for the conserved 

c 

electric current of a quark, which moves along a certain closed contour C. Then, performing 
the dualization of the so-extended action and summing up over monopole currents according to 
Eq. (||), we arrive at Eq. @ with F^ u replaced by + F*. Here, stands for the field 
strength tensor generated by quarks according to the equation d^F* = H- A solution to this 

equation reads F* = —gY? , where H e (x) = / da fiv (x(^))5(x — x(£)) is the so-called vorticity 

s c 

tensor current defined at an arbitrary surface S e (which is just the world sheet of an open dual 
Nielsen- Olesen string), bounded by the contour C. 

From now on, we shall be interested in the London limit, A — > oo, of the theories (||) and (|ll), 
where they admit an exact string representation. In that limit, the partition function of the 
theory (§) with external quarks reads 



Z 



VB^Ve smg -V6 TCS - exp 



- / d 4 x i (F, v 



— (8^9 — 2g m B fJ 



(12) 



In Eq. (0), we have performed a decomposition of the phase of the dual Higgs field 6 = smg - -|-# reg - , 
where the multivalued field 9 Bms '(x) describes a certain configuration of the dual strings and obeys 
the equation |22j 



e^x P d x d p e si ^ix) 



27r£ M1/ (x) 



(13) 



Here, £ Mi , stands for the vorticity tensor current, defined at the world sheet E of a closed dual 
string, parametrized by the vector x^). On the other hand, the field 9 ieg '(x) describes simply a 
singlevalued fluctuation around the above mentioned string configuration. 

The string representation of the theory ( O ) can be derived similarly to Ref . [25] , where this 
has been done for a model with a global Z7(l)-symmetry. One gets 



Z = VxJQVhwexpl - / d 4 x 



12r)' 



tt2 _,_ 2 1,2 



fill ^^pbV^y.V 



(14) 



where = — E MV , and H^x = d^h u \ + d\h pv + d u hx^ is the field strength tensor of a massive 
antisymmetric tensor field h pi/ (the so-called Kalb-Ramond field J53|). This antisymmetric spin-1 
tensor field emerged via some constraints from the integration over 6 re& ' and represents the massive 
dual vector boson. As far as the integration over the world sheets of closed strings, Dx^), is 
concerned, it appeared from the integration over 6* smg ' by virtue of Eq. (]T3|), owing to which there 
exists a one-to-one correspondence between # sm s- and £ M (£)- Physically this correspondence stems 
from the fact that the singularity of the phase of the dual Higgs field just takes place at the 
string world sheets. (Notice that since in what follows we shall be interested in effective actions 
rather than the integration measures, the Jacobian emerging during the change of the integration 



6 



variables (9 smg - — > which has been evaluated in Ref. will not be discussed below and is 

assumed to be included into the measure Vx^).) 

Finally, the Gaussian integration over the field hp U in Eq. (|14D leads to the following expression 
for the partition function ([13 



Z 



exp 



j dXf, j dy^D^\x - y) 
c c 



£>x /t (O ex P -(nrj) 2 J d 4 x J d 4 yt lfMU (x)D < £ ) (x - y)f,p U {y) 



(15) 



Here, D^(x) = ^,, Ki(m\x\) is the propagator of the dual vector boson, whose mass m, gen- 
erated by the Higgs mechanism, is equal to 2g m r], and K u 's henceforth stand for the modified 

IT 



Bessel functions. The details of derivation of Eqs. (fl^) and ( |15D can be found e.g. in Ref. [13|. 
Besides that review, the obtained string representation flip] ) has been discussed in various contexts 
12 , 17| , |24] , 25fl . Clearly, the first exponential factor on the R.H.S. of Eq. (|15|) is the 



in Refs. [10 



standard result, which can be obtained without accounting for the dual Nielsen- Olesen strings. 
Contrary to that, the integral over string world sheets on the R.H.S. of this equation stems just 
from the contribution of strings to the partition function and is the essence of the string repre- 
sentation. The respective string effective action describes both the interaction of the closed world 
sheets E's with the open world sheets S e 's and self-interactions of these objects. 

We are now in the position to discuss the bilocal correlator of electric field strengths in 
the model (|T2|) . Indeed, owing to the Stokes theorem, such an extended partition func- 
tion (which is actually nothing else, but the Wilson loop of a test quark) can be written as 



exp (-ffd A x^J, u 



where (. . .) f 



and (. . .) . m given by Eqs. (||) and @, respectively, 
expansion, we have in the bilocal approximation: 



JVa^exp (-Ses. a^f^ 



With S'eff. 



Applying to this expression the cumulant 



Z ~ exp 



x / d^l u (x)T,\ p {y){{f, u {x)h p (y))) a ^ 



(16) 



Following the SVM, let us parametrize the bilocal cumulant by the two Lorentz structures similarly 
to Eq. ®: 



((f^( x )hp(0))) a ^ = - 5 w 5 uX YD (x 



1 

+ 2 



8^ yx x 6 

Owing to the Stokes theorem, Eq. ( |TTD yields 



+ 



V 1 (x 2 ) . 



(17) 



Z * exp [-- J d*x j d A y [2^(x)E^(y)X> ((x - yf) + ((x - yf 

where 



(18) 



G x 



d\Vi(\). 



(19) 
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On the other hand, Eq. should coincide with Eq. flT5| ) divided by Z 



o 



just the standard normalization condition, encoded in the integration measures), i.e. it reads 



(which is 



Z = expl [ -Jd i xJ d A yD^\x - y) 



{A^l^lM + -jl(x)f u (y) 



x < exp 



8(7r V ) 2 J d*x J d 4 yD^( X -y)E; u (x)^(y) 



(20) 



where 



/X)x M (0(...)exp 


-(7r V ) 2 Jd 4 xJd 4 y^ u (x)D%\x-y)ZM 


JVx^exp 









(21) 



As it has already been discussed in the Introduction, in the case of zero temperature under study, 
dual Nielsen-Olesen strings, one should average over in Eq. (pip, form virtual bound states of 
vortex loops. The typical areas |E|'s of those are very small, and in the leading approximation can 
be disregarded w.r.t. the area |S e | of the world sheet of the long open string, which confines a test 
quark. Owing to this, the exponential factor, which should be averaged over vortex loops on the 
R.H.S. of Eq. Q2TiD, can be disregarded w.r.t. the first exponential factor in this equation, as well. 
Then, the comparison of the latter one with Eq. ([18]) straightforwardly yields for the function T> 
the following expression 



V(x 2 ) 



m 3 Ki(m\x\) 

47T 2 \x\ 



[22) 



whereas for the function T>\ we get the equation G [x 2 ) = AD${x), which leads to: 



V, (x 2 ) 



rn 



2tt 2 x 2 



Ki(m\x\) m 



\x\ 



K (m\x\) + K 2 (m\x\) 



(23) 



We see that in the IR limit \x\ ^> — , the asymptotic behaviours of the coefficient functions (22) 
and fl2"3|) are given by 



V 



m 



-m\x\ 



4v/2vt2 ( m \ x \y. 



(24) 



and 



m 



-m\x\ 



(25) 



2^7^ ( m | x |) 2 

For bookkeeping purposes, let us also present the asymptotic behaviours of these functions in the 
opposite case, |x| <C — . Those read 



V 



and 



m 



4ir 2 x 2 



(26) 



Pi 



IT 2 \x\ 



(27) 
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One can now see that according to the lattice data |3|, |j, the asymptotic behaviours 



and (p5) are very similar to the IR ones of the nonperturbative parts of the functions D and 
D\, which parametrize the bilocal cumulant (|l|) in the case of QCD. In particular, both func- 
tions decrease exponentially, and the function D is much larger than the function T>i due to the 
preexponential power-like behaviour. We also see that the role of the correlation length of the 
vacuum, T g , is the SVM, i.e. the distance at which the functions D and Di decrease, is played 
in the model fll2|) by the inverse mass of the dual vector boson, — . Moreover, the UV asymp- 
totic behaviours ( |2lp and fl2"7p also parallel the results of the SVM in QCD to the lowest order 
of perturbation theory. Namely, at such distances the function D\ also behaves as due to the 
one-gluon-exchange contribution. As far as the function D is concerned, it vanishes to the leading 
order of perturbation theory. Although this is not the case in our model (whose UV features are 
far from those of the asymptotically free QCD), the Pi-asymptotics ([27]) is nevertheless really 
much larger than that of the "D-function, given by Eq. (Pop. 

Hence we see that within the approximation when the contribution of vortex loops to the 
partition function (J2IJ) is disregarded completely, the bilocal approximation to the SVM is an exact 
result in the theory ([12]), i.e. all the cumulants of the orders higher than the second one vanish. 
Higher cumulants naturally appear upon performing in Eq. (|20|) the average fl2~I ) over vortex loops. 



However, this average yields important modifications already on the level of the bilocal cumulant. 
Namely, as we shall see in the next Section, it modifies the classical expressions (|22"D and (|23|). 

2.2.2 S77(3)-case 

Let us now turn ourselves to the bilocal cumulant of electric field strength tensors in the London 
limit of AP SU (3)-gluodynamics, where the partition function ( |TTj) of this theory takes the form 



Z 




exp< 



d x 



4 ^ 



2 3 
Z a=l 



?n^aB^) 



(28) 



Similarly to the >S[/(2)-case, in the model under study there exist dual Nielsen-Olesen-type strings. 
Due to that, in Eq. (^8|) we have again decomposed the total phases of the dual Higgs fields into 
the multivalued and singlevalued parts, 9 a = 6^ mg - + 6^ eg \ Here, the multivalued parts 6^ ing ''s 
describe a given configuration of the dual strings of three types. They are related to the world 
sheets S a 's of these strings via the equations 



2ttS^(x) = 2tt J da^,(x a (0)6(x - x a (0), 



(29) 



where x a = a^(£) is a four- vector parametrizing the world sheet S a . 

An external quark of a certain colour c = R, B,G (red, blue, green, respectively) can be 
introduced into the theory under study by adding to the initial action (^) the interaction term 
xsL^j^, where the vectors of colour charges read 



Q 



(R) 



Q 



Q 



o, 



V2' 2V3J ' V 2 ' 2^3, 

These vectors are just the weights of the representation 3 of *SU(3). In another words, those are 
nothing else, but the charges of quarks w.r.t. the Cartan subgroup [f/(l)] 2 , i.e. for every c, the 
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components of are just the charges of a quark of the colour c w.r.t. the diagonal gluons A 3 ^ 
and A* 

Applying further the Stokes theorem and the cumulant expansion in the bilocal approximation, 
we get for the partition function of the theory (^) with an external quark of the colour c the 
following expression: 



Z c ~ exp 



■ / d'x I d^ v (x)l? Xp {y){UiMf{M 



(30) 



Here, i,j = 1,2 denote the [£/(l)] 2 -indices, referring to the Cartan generators (T 3 ,T 8 ), and the 
average is defined as (. . -} a ^^ = (/£>a M exp (-S eS . a M ,f^ ) (. . , where S cS , and (. . .)_.„, are 

given by Eqs. (|7|) and (|TI]), respectively. Upon the SVM- inspired parametrization of the bilocal 
cumulant, 



Wf^O))) = 5»i 15» X S UP - 8, p 8 vX \D (x z ) 



1 

+ 2 



we can write for Eq. (^0) the following expression: 



D 1 [x 



(31) 



exp 



1 

24 



(Fx I d A y 2^ 2 S?.„(x)S e 



(IV \ 



(IV 



(y)D ((x - yf) + j;(x)j;(y)G ((x - yf)} } . (32) 



Here, we have denoted by G the same function as (|1^), but with the replacement T>\ — *■ Di, and 
used the fact that for every c, (Q^) 2 = |. 

On the other hand, one can derive the string representation for the partition function Z c . 
Indeed, the dualization of the action (0) with the term zQ^ / d^xa^j^ added, leads to Eq. (||) 
with F^ replaced by F^+F^j. Here, F^J stands for the field strength tensor of a test quark of the 



(IV 



(IV L ^yu^^u (IV 

colour c, which obeys the equation c^F^ = Q^jf, and thus can be written as FjfJ = — #Q 
Next, the summation over the currents of monopole Cooper pairs in the sense of Eq. ( |10D yields 
Eq. ( |TTD with the same extension of F^. In the London limit under study, the string representation 



of this theory (see Refs. JTI], ^6] for details) reads 



x 



x exp j-vr 2 J d 4 x J d 4 yD$ B (x - y) 



^{x^M + —j;( x )j;(y) 



(33) 



Here, mg = VQgmV is the mass of the dual vector bosons, which they acquire due to the Higgs 
mechanism. We have also introduced the notation S* = — 2s^H e with the following numbers 

1, which obey the 



,(«) 

'3 



IB) 
b 2 



JP) 



o, 4 R) 



SB) 



JG) 
b 2 



~ b 2 



sl c) 's: .s., 

- I e g (c) 

the following expression for the partition function (cf. Eq. 



6 i 



-(G) 



relation Q 1 -^ = ^e a s^K Taking into account that for every c, (s^) = 2, we eventually arrive at 



Zc = exp {-Syr 2 J d 4 x J d*yD$ B (x - y) 



^E* (x)£* (y) + — - j«(x)j«(y) 



48tt 



10 



x < exp 



(2nrj) 2 s^J d 4 x J d*yZ; u (x)DW(x-y)Zl u (y) 



K(0 



(34) 



with the average over vortex loops having the form 



(...) 



\a=l 



cxp 



(n V y J d±x J A^(x)^ B (x - y)E^(y) 



JVxoJOSl £ )exp [-(t^) 2 /^^" - (y) 



a=l 



Comparing now Eq. ( |3"2"D with Eq. (|3~3p, we see that in the approximation of very small vortex 
loops, |S a | <C |E B |, the functions D and Di are given by Eqs. (p^ ) and (|23"D with the replacement 
m — > mg. Besides that, it is obvious that the bilocal cumulant ( 51 ) is nonvanishing only for the 
gluonic field strength tensors of the same kind, i.e. for i = j = 1 or % = j = 2. Hence, for these 
diagonal cumulants, the vacuum of the AP S't/(3)-gluodynamics in the London limit does exhibit 
a nontrivial correlation length T a = —. 



Electric Field Strength Correlators in the Gas of Vortex 
Loops 



3.1 577(2). 



case 



To study the properties of vortex loops in the above considered theories, there is clearly no necessity 
to introduce external quarks. The field strength correlators can be studied afterwards, i. e. already 
after the summation over the grand canonical ensemble of vortex loops. Thus, let us first consider 
the theory (|T2|) with F* = 0. Upon the derivation of the string representation of such a theory, 
we are then left with Eq. (|i~4|) , where S B = 0. To study the grand canonical ensemble of vortex 



loops, it is necessary to replace in Eq. (|14|) by the following expression: 



N 
i=l 



da lxv {x i {^))8{x - Xi(g)). 



(35) 



Here, £ G [0, 1] x [0, 1] is a 2D-coordinate, and rij's stand for winding numbers. In what follows, 
we shall restrict ourselves to the vortex loops possessing the minimal winding numbers, = ±1. 
That is because, analogously to the 3D-case 0, pTjl , the energy of a single vortex loop is known 
to be a quadratic function of its flux, owing to which the existence of two vortex loops of a unit 
flux is more energetically favourable than the existence of one vortex loop of the double flux. 

Then, taking into account that the gas of vortex loops is dilute, one can perform the summation 
over the grand canonical ensemble of these objects, which yields the following expression for the 



partition function [17| : 



Z 



Vh 



exp 



d x 



12?/ 



TJ 2 _l_ n 2 h 2 



1m fiu 



2( COS 



"A 1 " 



(36) 



Here \h 



fjtv I 



and A 



is an UV momentum cutoff with L and a denoting the 



characteristic distances between vortex loops and their typical sizes, respectively. Clearly in the 



11 



dilute gas approximation under study, a <C L and A 3> a -1 . Also in Eq. Q5E|), C e~ 5 ° stands for 
the fugacity (Boltzmann factor) of a single vortex loop, which has the dimension (mass) 4 , with So 
denoting the action of a single loop. 

Note that the value of So is approximately equal to era 2 , where we have estimated the area of 
the vortex loop as a 2 , and a stands for the string tension of the loop, i.e. its energy per unit area. 
This energy can be evaluated from the action standing in the arguments of the exponents on the 
R.H.S. of Eq. ( ^TCP by virtue of the results of Ref. f28| and reads 



a 




(37) 



Here we have in the standard way |27[] set f°r a characteristic small dimensionless quantity in the 
model under study the value ^7=, which is of the order of the ratio of m to the mass of the dual 

Higgs field. Moreover, it has been assumed that not only ^> 1, but also In (^7) 1 
the last equality on the R.H.S. of Eq. (|3"7|) 



\9m 



i.e. 



origin of this logarithmic divergency is analogous to that, which takes place in 3D [14 



is valid with the logarithmic accuracy. The physical 

11 and is 



based on the fact that at the world sheet of a vortex loop the condensate of the dual Higgs field 
is destroyed, and the dual vector boson remains massless. 



The square of the full mass of the field following from Eq. (|36|) reads M 2 = m 2 +m 2 D = Q 2 rf 
Here, m D = i . , 

2 



+ •$$ is the full magnetic 



^2— is the additional contribution, emerging due to the screening of magnetic charge 

of the dual vector boson in the gas of electric vortex loops, and Q 
charge of the dual vector boson. 

To study the correlation functions of vortex loops, it is convenient to represent the partition 
function (|36D directly as an integral over these objects. This can be done by virtue of the following 
equality, 



exp 



d x 



12r/ ; 



H lv\ + 9mh 2 



m [11/ 



J PS^expj- (nr/) 2 J d 4 x J d 4 yS^ u (x)D ( £ ) (x - y)S^ v (y) + in j d 4 xh IMU S flv 



(3* 



in whose derivation it has been taken into account that d^h^ 



0. Indeed, owing to the Hodge 
decomposition theorem, the Kalb-Ramond field can always be represented as follows: = 
dfj,(p u — dv^fi + s^vXpdxipp. Clearly, in the original expression for the partition function, 



Z — If Vh^ exp< - J d A x 



1 



127/ 



tt'2 1 2 7 2 



1m /iv 



(39) 



gas 



the field decouples not only from S^ s (due to the conservation of the latter one), but also from 
tpn- The (/? M -field thus yields only an inessential determinant factor, which is not of our interest. 
Therefore this field can be disregarded, which proves the above statement. 

Substituting now Eq. (J38|) into Eq. (p6|), we can integrate the field out. This yields the 
desired representation for the partition function: 



Z = J VS^exp {- [(tt77) 2 j d 4 x J d 4 yS, u (x)D^(x - y)S, v (y) + V[S^\ 



(40) 
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where the effective potential of vortex loops, V, reads 



V[S, 



d x 





^ nA2 \s 1 + 


jyrA 2 !^! In 


2C 1 Vl \ 




2C 



\ 




• (41) 



It is straightforward to prove that the correlation functions of S^s, calculated by virtue of the 
representation (flO), are nothing else, but the correlation functions of vortex loops in the gas. This 
can be seen in the following way. Let us integrate the /i^-field out of the initial expression fl39|) 
for the partition function of the gas of vortex loops. This yields (c/. Eq. 



Z 



exp 



M 2 f d 4 x [ d 4 yE^(x)D^ 



x 



(42) 



This equation is now perfect to involve S^'s and demonstrate that their correlation functions are 
indeed equal to those of vortex loops. After that, the /independence can be restored back, so 
that we shall eventually again arrive at Eq. (^) with the substitution (|38|). In order to involve 
Sftv's, let us rewrite Eq. ( f42"D as follows: 



Z 



oo 

N=l iv • 



x ( / VS^S (S^ u - E£? 



The average here reads 



)exp [-(n V ) 2 J d'x J d i yS, v {x)D^\x-y)S, u {y)\) ^ 



{^(0}f = i 



(43) 



(O) 



N 

* L =11 / *Vi T>Zi(£)p[zi] £ O 



(44) 



ni=±l 



In this formula, the vector y^ describes the position of the world sheet of the z-th vortex loop 0, 
whereas the vector describes its shape, i.e. Xi(£) = Vi + z i(Oi Hi = I ^^(O- We have 
also denoted by fx a certain rotation- and translation invariant measure of integration over shapes 
of the world sheets of vortex loops. Note that it was just this average, which in the dilute gas 
approximation led from Eq. ( j3~9|) to Eq. ( p6|) (see Ref. |HJ for details). 

From the (5-function in Eq. (43) it is now clearly seen that the correlation functions of S^s 
are indeed equal to those of vortex loops. One can further represent this 5-function as an integral 
over the Lagrange multiplier, whose role, as we shall see immediately below, is just played by the 
Kalb-Ramond field: 



Indeed, normalizing the measure T>h^ v by the condition f] 



(45) 



J Vh^VS^exp -{nr]) 2 J d 4 x J d 4 ySf, v {x)D { ^ (x - y)S flu (y) - in J d A xh fn/ S fJ 

\k? B lv\ +9?n h lv } = L 



3 For brevity, we omit the Lorentz index. 

4 Owing to Eq. (|38|), this condition can be rewritten simply as J Vh^ exp j— J d A x 
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we get 



Z 



VS pv Vh pu exp exp 



{^f J d A x J d 4 yS, u (x)D<£\x-y)S^y) + 



+iir 



d xh[j,i/S hi/ 



2( cos 



' v P v | 



This is just Eq. with the substitution fl38|), which completes our proof. 

The correlation functions of vortex loops can now be calculated in the approximation when 
the loop gas is sufficiently dilute, namely its density obeys the inequality \S pi/ \ <C -p. Within 
this approximation, the potential (|4l| ) becomes a simple quadratic functional of S^s, and the 
generating functional for the correlators of vortex loops takes a simple Gaussian form. It reads 



Z[J, 



1 



P V\ 



Z[0] 



vs, 



exp 



(n V ) 2 { d 4 x [ d 4 yS pv {x)D^ 4 \x-y)S, v {y) + 



+ J d 4 x 
where 



TT^A 4 

4C 



•2C + 



. <?2 i J C 

u pu ~ °lLV t ~ J llV 



exp 



d 4 x / d 4 yJ pv [x)Q[x - y)J pv {y) 



(46) 



c 



G(x) = ^-(d 2 -m 2 )D^(x) 



7T 2 A 4 



(47) 



Next, since <9 M £^f = 0, the 5-function in Eq. P3| ) requires that d^S^ = as well. The Hodge 
decomposition theorem then leads to the following representation for S^ u \ S^ v = £ IJ ,uX P dx'fi P - Owing 
to this fact and the same theorem, the coupling / d 4 xJ pv S pv will be nonvanishing only provided 



ix 2 



that J pv = e^xpdxlp. This coupling then reads 2 / d 4 xI fJ T fll/ (p l/ , where T^ u (x) = dffi - S pu d 
On the other hand, substituting the above representation for J^ u into the R.H.S. of Eq. (pEBJ), we 
have 



Z[J t 



pv\ 



exp 



-2 J d 4 x J d 4 yI p ,{x)I u {y)T llv {x)Q{x - y) 



Thus, varying Z[J^ U } twice w.r.t. J M and setting then 1^ = 0, we get 

T^(x)Tx P (y) (ip„(x)(p p (y)) = -T lxX (x)G(x - y). 

Due to the rotation- and translation invariance of space-time, it is natural to seek for (ip u (x)ip p (y)) 
in the form of the following Ansatz: 5 up g(x — y). This yields the equation d 2 g = Q, whose solution 
reads 

c 



9(?) = ~^ dx2 - m ") J d'yD^ix - y)D<$(y), 



where D^\x) = D®{ 



x 



= . 2 2 - The last integral can obviously be rewritten as 

m=0 47r x 



Jd 4 zD^ 4 \z)D^{z- 



X 



(4£ 



As we shall see below, it will be necessary to know the more general expression, namely that for 
the integral 
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d i zD^(z)D^(z-x). (49) 
Its calculation is outlined in the Appendix, and the result reads 

^(D£\x)-D<$(x)). (50) 



m, 



Note that according to Eq. (|4"9"|), Eq. (^D|) should be invariant w.r.t. the interchange m <-> M. By 
noting that during this interchange m? D changes its sign, one can see that this invariance really 
holds. 

Setting now in Eq. fl50"D m = 0, we get Q -^r (Dq(x) — D$ D (x)\ which yields for the desired 
integral ( f48|) the same result with the substitution mj) -> M. Thus, the final expression for the 
function g reads 



/. 

The desired correlator of S,,,/s has the form 



g(?) = j^m^i 92 - ™ 2 ) K'w - ^ 4) (^)) • ( 51 ) 



ttMA 2 ) 2 

fiv s he 

(S llu (x)S Xp (y)) = e^ ua pe Xpia d x a d v (ipp(x)tp a {y)) 



and therefore 

(S IMU (x)S Xp (0)} = -e^ a /3Sxpj^^g(x) = 

= (SxuSpp - S up 5px) G{x) + (Spxd p d u + bypd^dx - 8 w d x d v - 5xvdpd p ) g(x), (52) 

where it has been used that d 2 g(x) = Q{x). 

This result can now straightforwardly be used for the calculation of the contribution of vortex 
loops to the bilocal cumulant ([H]). Indeed, applying to the average on the R.H.S. of Eq. ( p0|) the 
cumulant expansion in the bilocal approximation, we get: 



Z ~ expj- j d'x J d*yD$(x - y) [(^E^E^y) + \f^x) 3 ^{y) 



+ 



+32( 7 rr ? ) 4 / d 4 xd 4 yd 4 zd\D%\x - z)D$(y - u)K u (x)ZUy) ({E tlu (z)E Xp (u))) x 



Comparing this expression with Eq. (|IED, we see that owing to Eq. (0), the additional contribution 
of vortex loops to the cumulant (|T7|) has the form 



A((fp U (x)fx P (y))} a ^ = (4W)7 d 4 z J d 4 uD$(x-z)D£\y-u)x 
x {(5p X 5„p - 8 w 8 vX ) G(z -u)+ [5,pd{dl + S^Vffi - 8^1 - S^d^ g(z - u)} . 



5 Clearly, this result can also be obtained directly by making use of the method presented in the Appendix, 
which was done in Ref. lEql. 
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Let us further compare this intermediate result with Eq. QT7D and take into account that 



(x - y\V x ((x - y) 2 ) = ((x - y) 2 ) 



where the function G is defined by Eq. (|T9|). This leads to the following system of equations, 
which determine the contributions of vortex loops to the functions D and G: 



u 



AV((x-y) 2 ) = (W)7 *z J d i uD^(x-z)D^(y-u)g(z- 

J d'z J d 4 uD%\x - z)D$(y - u)g(z - u) 



J d 4 uD£\y-u)D$(x-u). 



AG ((x - yf) = (87rg mV 
Inserting now Eq. into Eq. (|53|), we get 



(53) 
(54) 



AVUx-yf 



(% m r/ 2 ) 2 C 



A 4 



By virtue of the Appendix, this yields 



AV x 



m 
4^2 



M m 
— i^i(M|x|) - — ATxHxl 



x 



Adding this result to Eq. fl22|), we finally obtain for the finction P the following full result: 

m 2 MATi(M|x|) 



Analogously, inserting Eq. 



AG ((x - y) 2 ) = C 



pftlll 

U) into Eq. 

igmV \ 



X 



4tt 2 
we have 



x 



(55) 



v A 2 M 

or further by virtue of the Appendix, 

AG (x 2 



d\D$(y 



u 



( 4 ), 



(^)'^)-^£>(z). 



ttM\x\ 



Together with Eq. (p3"D, this yields the following full result for the function T>\\ 



(x 2 ) 



m 



D 



n 2 M 2 



+ 



m 



x\ 



2n 2 Mx 2 



K x {M\x\) M 



x 



(K (M\x\) + K 2 (M\x\)) 



(56) 



It is worth noting that the functions AT) and AT>\ contain the terms exactly equal to Eqs. fl22|) 
and (|23|), respectively, but with the opposite sign, which just cancel out in the full functions (|55|) 
and dHI). We also see that, as it should be, the functions (|55|) and (|56|) go over into Eqs. (|^) and 
(p3|), respectively, when mc — >• 0, z.e. when one neglects the effect of screening in the ensemble 
of vortex loops. An obvious important consequence of the obtained Eqs. (|55f) and (p6|) is that 
the correlation length of the vacuum, T g , becomes modified from — (according to Eqs. ( p2|) and 
(f23|)) to -j. (It is worth emphasizing once more that this effect is just due to the Debye screening 
of magnetic charge of the dual vector boson in the ensemble of electrically charged vortex loops, 
which makes this particle more heavy, namely enlarges its mass from m to M.) Indeed, it is 
straightforward to see that at |x| 3> ^g, 
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and 



^ 4V2ttI (M|x|)§ 



pfuii _^ rn 2 D (ml) 2 e~ M N 

1 _9. _ld ' 



tt 2 M 2 |x| 4 2^7rl (Af|i|)f 

It is also remarkable that the leading term of the IR asymptotics of the function £>J u11 is a pure 
power-like one, rather than that of the function T>i, given by Eq. fl25|). Another nontrivial result 
is that the screening does not change the UV asymptotic behaviours of the functions (|22|) and 
(f23|), i.e. the UV asymptotics of the functions (|55D and ( |56"D are given by Eqs. fl26|) and 
respectively. 

Finally, it is worth remarking that due to the modification of the D-function, one could expect 
the appearance of some change in the string tension of the open dual string world sheet S c . 
However, by virtue of the general formula expressing the string tension via the P-function [2S|, 



a = 4Tg J d 2 zT>{z 2 ) ) one can check that this is not the case, i.e. the string tension of S e is 
independent of whether we account for screening in the gas of vortex loops or not. The reason 
for that becomes clear from the resulting expression for a. It reads IQnrj 2 In - with c standing for 
a characteristic small dimensionless quantity, and thus depends only on rj, which is not affected 
by screening. Similarly to Eq. (p7|), setting for c the value we see that the string tension 
of S e is in the factor 16 larger than the string tension of a vortex loop. Clearly, that is due to 
the factor 4 standing in the linear combination of and in E M! , (cf. Eq. (|1"5|)). However, 
the coupling constant of the next-to-leading term in the gradient expansion of the nonlocal string 
effective action standing in the second exponential factor on the R.H.S. of Eq. ( pT5[ ) (the so-called 
rigidity term) does depend explicitly on the magnetic coupling constant and therefore changes due 
to the screening. Indeed, by virtue of the results of Refs. |2S|, [T3|] , one can see that for the same 



world sheet E e , this coupling constant without taking screening into account reads j^-y, whereas 
in the presence of screening it goes over to IS-, as it could be intuitively expected. 



3.2 5'^(3)-case 

In the present Subsection, we shall extend the above results concerning the electric field strength 
correlators in the gas of vortex loops to the case of AP S'[/(3)-gluodynamics (p8|) . When the 
dual Nielsen-Olesen strings in this theory are considered as noninteracting objects, the string 
representation of its partition function is given by Eq. (|33|), where one should set = 0. 
Integrating out the coordinates of one of the three world sheets (for concreteness, a^(0)> we can 
write the so-obtained expression for the partition function as 

Z = JvxKOVxKOx 



xexp | -2M 2 / d*x I d'y K V {*KM + K^lM + K»( X KM ~ v)\ ■ (57) 



In order to proceed from the individual strings to the grand canonical ensemble of vortex loops, 
one should replace (x), a = 1,2, in Eq. (pTF]) by 
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N 

i=i J 



Here, n°'s stand for winding numbers, which we shall again set to be equal ±1. Performing 
such a replacement, one can see the crucial difference of the grand canonical ensemble of vortex 
loops in the model under study from that in the AP SU (2)-gluodynamics, studied in the previous 
Subsection. Namely, the system has now the form of two interacting gases consisting of the vortex 
loops of two kinds, whereas in the SU (2)-case the gas was built out of vortex loops of the only 
one kind. 

Analogously to that case, we shall treat such a grand canonical ensemble of vortex loops in 
the dilute gas approximation. According to it, characteristic sizes of loops are much smaller than 
characteristic distances between them, which in particular means that the vortex loops are short 
living (virtual) objects. Then the summation over this grand canonical ensemble can be most 
easily performed by inserting the unity 



i = / vs;j (s; v - E«r) (5* 



into the R.H.S. of Eq. ( p7|) (where E" is replaced by E°^ as ) and representing the (^-function as 



an integral over the Lagrange multiplier (cf. Eqs. (|43|) and fl45|)). Then, the contribution of N 



vortex loops of each kind to the grand canonical ensemble takes the following form: 



Z 



y-ia gas 
^/xv 



vs; v v\%x 



x exp 



2{nr,f I d'x I d"y [^(x)^(y) + tW^M + 3U*)<Sj,(l/)l 4S> " v) 



hi/ \ 



/IV 



/IV \ 



/IV * 



/IV * 



-i J d* x \% (s; u - s» 



gas 

/IV 



(59) 



After that, the desired summation over the ensemble of loops is straightforward, since it technically 
parallels the one of AP SU (2)-gluodynamics. We have 



n 

0=1 



00 C N I ft 
1 + E (^P (i J d A x\ a ^l 



exp | 2C J d% 



cos 



KM 



A 2 




KM 



A 2 



(60) 



where for every a, the average (. . -}{ x o.^n is given by Eq. (fH|). Here, it has been naturally 

assumed that the vortex loops of different kinds have the same fugacity £ oc e~ 5 °, since different 
flsmg.jg en ^ er initial partition function ( p8|) in the symmetric way. Clearly, the action So of 
a single loop can be estimated analogously to how it has been done in the previous Subsection 
for the SU(2)-case. In Eq. (|B0|), we have also introduced a new UV momentum cutoff A = 
(^> a -1 ), where a again denotes a typical size of the vortex loop, whereas L stands for a typical 
distance between loops, so that in the dilute gas approximation under study a <C L. Finally in 

Si 



Eq. (|60D, we have denoted 



A /IV 



/IV 
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Next, it is possible to integrate out the Lagrange multipliers by solving the saddle-point equa- 
tion following from Eqs. ([59]) and (|60D , 



A fJ,U 



sin 




iA 2 

■ 2 £ <V- 



After that, we arrive at the following representation for the partition function of the grand canon- 
ical ensemble of vortex loops: 



Z = J VS; u exp {- [2(n V fJ d'x J d*y [S^x) S^y) + ^{x)^{y) + S^S^y) 



x 



D$(x-y) + V 



ci l 

7T 



V 



C2 
7T 



(61) 



which owing to Eq. flog ) is natural to be referred to as the representation in terms of the vortex 
loops. In Eq. (|6l|) , the effective potential of vortex loops is given by Eq. fldip . 

Next, to get the Debye masses, corresponding to the two interacting gases of vortex loops, it 
is necessary to get the respective sine-Gordon theory of the two Kalb-Ramond fields. In order to 
derive it, let us first diagonalize the quadratic form in square brackets on the R.H.S. of Eq 
which can be done upon the introduction of the new integration variables <S^„ = ^ {^Ju/ 

\ ( s lv - s %)- After that > E q s - flU) and © y ield 



fJ,U 



and S 2 U 



Z = J VS^VX^^-2^) 2 J d'x J d 4 yS; u (x)D^ B (x - y)S%{y)- 



+2C/ 



d 4 x 



cos 



A 2 




cos 



A 2 



(62) 



Here, we have introduced the two Kalb-Ramond fields as the following linear combinations of the 
Lagrange multipliers: = 4= + A 2 ^) and h 2 ^ = X 1 ^ — X 2 LU . The partition function of the 
desired sine-Gordon theory can then be obtained from Eq. (|62"D by making use of the following 
equality: 



J VS^expL W^f J d 4 x J d i yS%{x)D^ B {x-y)S; v {y)+i j ^xh%^ 



exp 



2tt 2 



d 4 x 



12r) 



(63) 



Similarly to the analogous equality fl38|), the equality (|63| ) can easily be proved by noting that due 
to the Hodge decomposition theorem and the equation d^S^ = (which follows from Eq. (]58[) 
and conservation of Sf„f as ), = 0. Substituting further Eq. ( |6"3"D into Eq. (p2] ) and performing 



h" 



/IV 



the rescaling -^j= — > h^, we arrive at the following representation for the partition function of the 
grand canonical ensemble of vortex loops in terms of the local sine-Gordon theory, equivalent to 
the nonlocal theory (|6l|): 
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Z = J Vh^ v exp | - / d 4 x [ JL (H; uX ) 2 + 3 -gl (h%) ' - 



-2C 



7T 



COS 



\A 2 \/2 



\/3/j 



+ cos ( = V^/iL — h 2 nj 

\A 2 V2 ^ 



(64) 



The full masses of the Kalb-Ramond fields can now be read off from Eq. ( |6"4]) by expanding the 
cosines up to the quadratic terms. The result reads M 2 = m 2 B + m 2 = Q 2 r] 2 , where mi = ^?a/3C, 
^2 = ^r?v^C are the Debye masses, and we have introduced the full magnetic charges Q\ = 



Next, Eq. (|61| ) can be used for the evaluation of correlators of vortex loops, which due to 
Eq. (|58| ) are nothing else but the correlators of S^'s. Those are again calculable in the approx- 
imation of a dilute gas of vortex loops, 5" <C Within this approximation, the generating 
functional for correlators of S^s reads 



ja 



PS^ exp 



2{-Krjf / d 4 x / d i ySUx)D^(x-y)SUy) + 



I 



+ d 4 x 



A* 

2C V3 



liu " [IV 



where J" is a source of S®, and J* = J* ± J 2 Apart from an inessential constant factor 
(which can as usual be referred to the integration measure and eventually drops out during the 
calculation of the correlation functions), we thus get for the generating functional the following 
expression: 



Z 



ja 



exp 



d 4 x J d 4 y [Ji,(x)J^,(y)g 1 (x - y) + J pv (x)J pv {y)Q 2 {x - y)] \ , 



where Q a (x) = (d 2 — m 2 B ) D^ a (x). Owing to the conservation of 5" 's and the Hodge decom- 
position theorem, we again have the following representation for 5" : S° = e pu \ p d\Lp a p . Therefore, 
due to the same theorem, J" = e pu \ p d\I p , which yields 



Z 



ja 



exp< —2 / d 4 x / d 4 y 



i;{x)Il(y)T^{x) (&(x -y) + G 2 (x - y)) + 



+2ll(x)I 2 (y)T tlu (x) (&(z -y)- Q 2 {x - y)) 

On the other hand, the coupling / d 4 xJ pv S pv can be written as 2 / d^xI^T^ip®. Thus, varying 
twice w.r.t. J"'s and setting then 1° = 0, we arrive at the following system of equations: 



Z 



ja 



T pi/ (x)T Xp (y) (<pl(x)<p (y)) = T pi/ (x)T Xp (y) (<p u (x)<p (y)) = -T M (x) (Gi(x - y) + G 2 (x - y)) 



T pu (x)T Xp {y) ((pl(x)tp 2 (y)) = -T pX (x) (g x (x - y) - Q 2 {x - y)) 
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Adopting for the correlators of ip^s the following Ansatze, 



we get: 



c 



2A 4 



rn 



B 



± - D<f\x)) ± ^ «(*) - 

This result makes the choice of notations u f±(x) n quite natural. Finally, the desired correlators 
of vortex loops read 



+ (Sn\d p d u + S up d p d x - S pp d x d u - S^d^dp) f+(x), 



(65) 



+ (Sp X dpd u + S up d p d x - S pp d x d u - b vX d p dp) /_ (x) . (66) 
This result can now be applied to the calculation of the contribution to the correlator (j31~l), 



brought about by the vortex loops. Indeed, applying to the average over vortex loops, standing 
on the R.H.S. of Eq. (|34]) the cumulant expansion in the bilocal approximation, we have due to 
Eq. ©: 



A((^(x)^(y) ; 

= -247rW4 c) 4 c) / d'z J d%D%l(x - z)D$ B (y - u) {S%{z)S\ p {u)] 
Taking further into account the equalities 

'sl v {x)sl P {y)) = {slMsl P {y)) = {slMslM] 

and the facts that for every c, (s^) = 2, s^f 1 s^ 1 + sf^ + s^s^ = —1, we can write 

*M e) (S a P v(z)S b Xp (u)) = 2 W) ~ (fyWW)) ■ 

This leads to the following system of equations: 

AD ((x - yf) = 487T 2 ^ 4 / *z J d 4 uD$ B (x - z)D$ B (y - u)G 2 (z - u), 
AG ((x - yf) = 96ttW/ d 4 z J #uD$ B {x - z)D<$ B {y - u) (f + (z - u) - f_(z - «)) 



where G is given by Eq. (|19D with the replacement T>\ — ► Dx. Carrying now the integrals out 
analogously to how it was done in the previous Subsection for the SU (2)-case, we get 
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AD (x 2 ) = m% (D 



(4) 
M 2 



X 



£)(4) f 



X 



AG x 



£>(4) / 



X 



— )D i4) (x) 



Together with the old expressions for the functions D and Di (given by Eqs. (H) and (H), 
respectively, with m —>■ tub), which did not account for the screening effect in the gas of vortex 
loops (], we finally obtain that D fu11 and Df 111 are given by Eqs. ( |55|) and (|56"1), respectively, with the 
replacements m — > tub, mu — * rri2, and M — > Mi. Therefore the whole discussion, following after 
Eq. (|56D, remains the same modulo these replacements. In particular, when ni2 vanishes, i.e. one 
disregards the effect of screening, the old expressions for the functions D and D\ are recovered. 



4 Modifications of the Propagators of the Dual Vector 
Bosons due to the Screening in the Gas of Vortex Loops 

In the present Section, we shall investigate the influence of screening of the dual vector bosons 
in the gas of vortex loops to the propagators of these bosons themselves. Let us start with the 
SU{2)-case ( |T2"D by studying the Wilson loop of a test magnetic particle, whose charge is in the 
factor n larger than that of the dual Higgs field. Such a Wilson loop has the form 



(W(C)) = ^exp (2ig m n J d^xB^jX) ^ exp 



2g m n d x Id yj^x) (B fl (x)B u (y)) j v (y) 



(67) 

where j^x) = § <ix Ai (s)5(a; — x(s)), and the average is defined by the partition function (|T2|) with 

c 

= 0. Clearly, in the derivation of the last equality on the R.H.S. of Eq. ( |6"TD we have used the 
cumulant expansion in the bilocal approximation. On the other hand, one can derive the string 
representation for this Wilson loop, and the result reads f24j| : 



(W{C)) = exp 



2g 2 m n 2 f d*x f d*yj,(x)DW(x-y)j,(y) 



X Z / Vx ^ ex V -( nr l) 2 J <fx J d 4 yi:^(x)D^\x - y)S^(y) + i J d i xJ fll/ E fll/ 



where 



J^u(x) = nne^xp I (Pyjx(y)dP \D$(x - y) - D^\x - y) 



Using again the bilocal approximation for the average over vortex loops, we have 



(W(C)) = exp 



-2g 2 m n 2 J d 4 x J d A yj,{x)D$(x-y)j»{y)- 



(6£ 



6 It is remarkable that these expressions again become exactly cancelled by the corresponding terms in AD and 
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d 4 x / d 4 y,J^ u (x) (S^ u (x)S Xp (y)) J\ P {y) 



(69) 



where the correlation function of the vortex loops, (Sfj, u (x)Sx p (y)), is given by Eq. (|52|). Inserting 
now Eqs. (|52|) and (^) into Eq. (^), comparing the latter one with Eq. (|6"TD, using the conservation 
of jp and the fact that d 2 g = Q, we get after straightforward calculations: 



{B^x)B v (0)) = 6, 



+ 9 - f *zj d\ (DjJ(i -z)- D { *\x - z)) Q{z - u)d 2 (D$(u) - nfiu)) 



(70) 



Substituting now here the explicit expression (|47f) for the function Q and using the Appendix for 
the calculation of the resulting integrals, we eventually arrive at the following simple expression 
for the propagator of the dual vector boson: 



1 



(B,(x)B u (0)) = S^jp (rnlD^tx) + m 2 D 



2n(4), 



M 



X 



(71) 



This equation replaces the classical expression 5 llll D < ^{x), one has without taking screening into 
account. Clearly, as it should be, this old expression is reproduced upon taking the limit mo — > 
in Eq. ([71]). 

Let us now consider the SU(3)-case fl28D , where the magnetic Wilson loop has the form 



(W(C U C 2 ,C 3 )) = ^exp (2ig m n J d 4 xB / J / 



exp 



-2g 2 m n 2 J d'x J Jyjftx) (B;(x)B b u (y)) j h v {y) 



(72) 



Similarly to Eq. (|9]) , we have used here the fact that the monopole charges are distributed over the 
root lattice. Therefore, we have set j M = e a j®, where j^(x) = § dx^(s)5 (x — x^ a \s)^ with the 



B„e„ 



contour C a parametrized by the vector x^'(s). We have also introduced the notation B^ _ 
and used the cumulant expansion in the bilocal approximation. In what follows, we shall be 
interested in the expression for the propagator (b^x) B b u {Qi)^ . To derive it, let us again consider 
the string representation for the Wilson loop under study. It turns out to be analogous to the one 
we had in the <SC/(2)-case and reads: 



(W(C 1; C7 2 ,C3)) = exp 



-3g 2 m n 2 J d'x J d'yj^D^ix - y)j;(y) 



x 



x 



i J Px^O^E^exp [-(n V ) 2 J d'x J d'yE^D^x - y)S^(y) + i J d±xJ^ v 



where 



J a (x) = irne^ I d A yf x (y)cf p D<£> B {x - y) - D^(x - y) 



(73) 



Applying the cumulant expansion, we have in the bilocal approximation: 



23 



{W{C h C 2 ,C 3 ))=exp 



-3g 2 m n 2 / d*x / dSfJx)D^(x-y)fJy) 



d*x / d*yj; u {x)(S« y {x)S b Xp (y))J b Xp (y) 



(74) 



where the correlation functions of vortex loops, (S^ u (x) S\ p {y)j , are defined by Eqs. ( |G5] ) and 
(|B(f). These equations, once being inserted into Eq. ( [74] ) together with Eq. ([T3|), yield upon some 
calculations and comparison of the result with Eq. (|72|): 



BUx)Bl(0)) = (Bl(x)Bl(O) 




( m B 



VmiMi 



D (4) 



f_mB_ 

\m 2 M 2 



D (4) 



M x 2 + M| 
MlMl 



D 



(4), 



(75) 



mi Mi 



Vm 2 M 2 



2m^ 



2 \r2 D ) ( X ) 



M(M 2 



■ (76) 



These equations represent the desired modifications of the propagators of the dual vector bosons 
in the model (^) due to the Debye screening of the magnetic charges of these bosons in the 
gas of electric vortex loops. As one can see, in the limit when this effect is disregarded, i.e. 
mi,m 2 tub, the diagonal propagators (|75| ) go over to the classical expression 5^ U D^ 



. x 



whereas the off-diagonal one, (|76|), vanishes. This means that the (quantum) effect of screening 
leads in particular to the appearance of the nontrivial correlations between dual vector bosons of 
different types. 



5 Conclusion 

In the present paper, we have investigated field correlators in the AP SU (2)- and SU (3)-theories. 
These correlators are of the two types. First of them are the correlators of electric field strengths, 
which are the most important ones. That is because they correspond to the gauge-invariant 
correlators in the real non-Abelian theories, which play the major role in the SVM and owing 
to that are widely used in the phenomenological applications. In AP theories, such correlators 
have up to now been evaluated only classically. In the present paper, we have improved on these 
calculations by evaluating the contributions to the correlators, brought about by the screening of 
the dual vector bosons in the gas of virtual vortex loops. This effect is essentially quantum as 
well as such a gas itself. In this way, it has been found that the correlation length of the vacuum 
in the models under study becomes modified from the inverse mass of the dual vector bosons, 
those acquire by virtue of the Higgs mechanism, to their inverse full mass, which takes also into 
account the effect of Debye screening. Besides that, in one of the two coefficient functions, which 
parametrize the bilocal correlator of electric field strengths within the SVM, there appears also 
a nontrivial power-like IR part, which was absent on the classical level. It has also been checked 
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that in the limit when the effect of screening is disregarded, the obtained novel expressions for the 
bilocal correlators in the AP theories go over to the classical ones, as it should be. It has also been 
discussed that the found modifications of the bilocal correlators do not affect the string tension, 
since this quantity depends only on the v.e.v. of the dual Higgs field and not on the mass of the 
dual vector bosons. Contrary to that, the coupling constant of the so-called rigidity term changes 
due to the screening. 

The correlators of the second kind, we have discussed in this paper, were the propagators of the 
dual vector bosons themselves. It turned out that (contrary to what happened to the correlation 
length of the vacuum in electric correlators) in the magnetic propagators the effect of screening 
does not lead simply to the change of the mass from the pure Higgs to the full one. There rather 
appear the expressions quite of a novel form, which, however, also go over to the classical ones 
when the effect of screening is disregarded. Besides that it has been found that in the S77(3)-case, 
screening leads also to the appearance of the nonvanishing correlations between the fields of the 
dual vector bosons of different kinds. This effect is a purely quantum one and disappears in the 
classical limit, when the screening is disregarded. 

In conclusion, the obtained results shed some light to the vacuum structure of the AP theories 
and give a new field-theoretical status to the SVM of QCD. 
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7 Appendix. Calculation of the Integral Qz£9l ) 

In this Appendix, we shall present some details of calculation of the integral fl49|) . Firstly, owin 
to the definition of the functions D$ and Dj$, we have: 



J*zD<$(z)D®(z-x) 



d*p r d 4 q f , 4 e ipz e^* - *) 
a z- 



(27r) 4 7 (2tt) 4 ./ p 2 + m 2 q 2 + M 2 
d A p e tpx 



(2TT) 4 (p 2 + m 2 )(p 2 + M 2 ) 
Next, this expression can be rewritten as 



,4 +°° +5° -l +°° +°° -am 2 ~/3M 2 - Al x l m 

P-J da\ d 0e^«<^ a HV+^ = J— da df— — — . (Al) 

2tt 4 J J H (An) 2 J J H (a + 3) 2 K ' 

v ' v ' v rj 

It is further convenient to introduce new integration variables a G [0, +oo) and t G [0, 1] according 
to the formulae a = at and j3 = a(l — t). Then, the integration over t yields for Eq. (A.l) the 
following expression: 
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+ 0O 

o 



(47TOTd) 

Such an integral can be carried out by virtue of the formula 

+00 /a\ I 

y* z"- 1 e---' ja! da: = 2 (2^7) , 3fy3 > 0, K 7 > 

and the result has the form of Eq. ( 150 ) from the main text. 



0. 
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